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Abstract 

We review some work done by us calculating matrix elements for Kaons. 
Emphasis is put on the matrix elements which are relevant to predict non- 
leptonic Kaon CP violating observables. In particular, we recall our results 
for the Bk parameter which governs the K° — K° mixing and update 
our results for e' K including estimated all-higher-order CHPT corrections 
and the new results from recent analytical calculations of the AI = 3/2 
component. Some comments on future prospects on calculating matrix 
elements for Kaons are also added. 
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1 Outline 



Here we review some work done by us in [1, 2, 3, 4, 5] devoted to the calcula- 
tion of matrix elements of Kaons. The main motivation for these calculations 
is the reduction of the hadronic uncertainty in those matrix elements which are 
necessary to test the Standard Model and unveil beyond the Standard Model 
physics. 

We will limit our discussion to matrix elements related to CP-violating effects 
involving Kaons: namely, the -B^-parameter governing the indirect CP-violation 
in the K° — K° mixing [1, 4] and the direct CP- violating parameter e' K [3, 5]. We 
update e' K by including the non-Final State Interaction (non-FSI) contributions, 
the new full isospin breaking corrections [6] as well as new results for some of the 
inputs used. We also use the results from the analytical results for the AI = 3/2 
component [5] to substitute the ones used in [3], these two are, incidentally, in 
numerical agreement. 

We don't discuss the light-by-light hadronic contributions to the muon g— 2 [7] 
nor the matrix elements of Kaons needed to quantify the electromagnetic mass 
difference [8]. Both calculations were done using techniques similar to those 
discussed here and can be improved along the lines discussed in the last section. 

2 Motivation and Notation 

We study here matrix elements related to two aspects of CP violation in Kaon 
physics; namely, 

a) The Bk parameter, defined as [C(u) is a Wilson coefficient] 

C(u) (K°\ J d 4 x Q AS =2(x)\K°) = ^B K f K m 2 K ; 

Qas= 2 (x) = 4L fl (x)L^x) and (1) 
2L^x) = [sY(l-l 5 )d](x). 

which enters in the indirect CP violating parameter Ek- The B^ parameter is 
an important input for the analysis of one of the Cabibbo-Kobayashi-Maskawa 
(CKM) unitarity triangles -see [9] for more information- 

and b) The direct CP violation in K — * im decays parameterized through e' K , 
defined as 

y/2e' K _ A [Kl — > (tttt) J=2 ] A [K s -> (tttt) J=2 ] 
e k A[K L ^(irn) I=0 ] A[K s ^(mr) I=0 }' 

In the isospin symmetry limit, K — > im invariant amplitudes can be decomposed 
into definite isospin amplitudes as [A = —iT] 

tA[K°^n\°] = ^ e **»-?^ e *, 
L V3 



1 



iA[K° 



\/3 VQ 



(2) 



with 5 and 5 2 the FSI phases and under very reasonable approximations, one 
can get 



£ K 



ie i(52-S ) Refl2 

a/2 Re ao 



Im a Im a 2 



Re ao Re a 2 



(3) 



To lowest order (LO) in Chiral Perturbation Theory (CHPT), i.e. order e°p 2 
and e 2 p°, strong and electromagnetic interactions between n, K, r\ and external 
sources are described by 



£(2) = _JL tr (uX + X+ ) + e 2 C 2 tr (QUC^) 



(4) 



with U = u 2 = exp(A a 7r a /F ) and = iv) (D^U)^ . \ a are the Gell-Mann matri- 
ces and the n a are the pseudo-scalar-mesons n, K, and rj. Q — diag(2/3, —1/3, —1/3) 
is the light-quark-charge matrix and x+ — u^x u ^+ u X^ u an d X — 2-Bo diag(m u , rrid, m s ) 
collects the light-quark masses. To this order, f n = Fq = 87 MeV is the pion 
decay coupling constant. Introductions to CHPT can be found in [10]. 

To the same order in CHPT, the chiral Lagrangians describing |AS| = 2 and 
| AS' | = 1 transitions are respectively 



,(2) 

'|AS|=2 



CF* G AS=2 tr (A 32 u^) tr (A 32 u M ) + h.c. 



(5) 



and 



,(2) 

"|AS|=1 

+CF* 



= CF 6 e 2 G s tr(A 32 utQu) 
G , 8 tr(A 32 uX) + G / 8 tr(A 32 x + ) 



+G 27 ^ fci tr(A ij U M )tr(A kl u^) 



+ h.c. 



(6) 



Here A {j = uXijvt and (Xij) a b = 5 ia 5 jb . 

3 Gp 



C = 



5 V2 VudKs 



-1.07 x 10~ 6 GeV 



-2 



(7) 



The normalization C is a known function of the PF-boson, top and charm quark 
masses and of the CKM matrix elements. The SU(3) x SU(3) tensor t^' kl can be 
found in [11]. 

In the Standard Model (SM), ImG 27 vanishes and ImG 8 and lmG E are pro- 
portional to CP- violating phases -to Imr with r = —X t /X u and Aj = V ic iV* s where 
Vij are CKM matrix elements. 
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At this order one gets 

B X K = \ Gas=2 , (8) 

'Reao\ (2) _ ^ 9Re G 8 + G 27 
:a 2 ) ~ IOG27 

/Imaoy 2 )^ ImGs (10) 
\Rea J ReG s + G 27 /9' V ; 



/Reao\ (2) ^ 9ReG 8 + G 27 , , 

UeaJ 10G 27 ' W 



and 



/Ima 2 y 2 )^ 3 F 2 Im (e 2 G g ) _ 
\Rea2/ 5m 2 K -ml G 27 

where we have disregarded some tiny electroweak corrections proportional to 
Re(e 2 G £ ). 

The large number of colors (N c ) predictions for the couplings in (5) and (6) 
come from factorisable diagrams, one gets 

G A s=2 = G 8 = G 27 = 1; and G E = . (12) 

Our goal in [1, 2, 3, 4, 5] was to calculate the next-to- leading order (NLO) in 
l/N c contributions to these couplings as well as the chiral corrections to Bk- 
Due to the lack of space, we will not go into details about the technique used 
to calculate hadronic matrix elements in those references. All the details on the 
X-boson method and on the short- distance matching were given there. Just to 
comment here that, in general, we compute a two-point function 

n,,(g 2 ) = (13) 
% J tfxe iq - x (0|r(i^ t (0)P J -(a:)e irLD ) |0) 

in the presence of the long-distance effective action of the Standard Model Tld- 
The pseudo-scalar sources Pi(x) have the appropriate quantum numbers to de- 
scribe K — > 7r transitions. The effective action Tld reproduces the physics of the 
SM at low energies by the exchange of colorless heavy X-bosons. To obtain it 
we make a short-distance matching analytically, which takes into account exactly 
the short-distance scale and scheme dependence. We are left with the couplings 
of the X-boson long-distance effective action completely fixed in terms of the 
Standard Model ones. This action is regularized with a four- dimensional cut-off, 
He- Our X-boson effective action has the technical advantage to separate the 
short- distance of the two-quark currents or densities from the purely four-quark 
short- distance which is always only logarithmically divergent and regularized by 
the X boson mass in our approach. The cut-off fic only appears in the short- 
distance of the two-quark currents or densities and can be thus taken into account 
exactly. 
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Taylor expanding the two-point function (13) in q 2 and quark masses one can 
extract the CHPT couplings Gas=2, Cg, G27, ■ • • and make the predictions of the 
physical quantities at lowest order. One can also go further and extract the NLO 
CHPT weak counterterms needed for instance in the isospin breaking corrections 
or in the rest of NLO CHPT corrections. 

After following the procedure sketched above we are able to write B K) ReG 8 , 
ImG 8 , and G 2 j as some known effective coupling [3, 4] \gi\ 2 (M x , fic, • • •) times 



where Uppab{Q 2 , <f) is a four-point function with AB being either L^L^ or 
L^R^; and W 1 are left and right currents, respectively, and Q = \px\ is the 
X-boson momentum in Euclidean space. The case of Im(e 2 G j E) can be written 

as [5] 



Therefore in this case one can use data as done in [5, 12, 13]. In [14], a Minimal 
Hadronic Approximation to large N c was used to saturate the relevant two-point 
function. 

We split the to 00 integration in (14) and (15) into a long-distance (LD) 
piece (from to /i) and a short-distance piece (SD) (from /i to 00). The SD piece 
we do using the Operator Product Expansion (OPE) in QCD and the result 
is thus model independent. For the long-distance piece we used data for Ulr 
and, as a first step, the ENJL model [15] for the four-point Green's functions. 
The good features and drawbacks of the ENJL model we used have been raised 
several times in [3, 7]. The most important drawback is that it does not contain 
all QCD constraints [16]. There is work in progress to substitute it by a ladder 
resummation inspired hadronic model that in addition to the good features of 
the ENJL model (contains some short- distance QCD constraints, CHPT up to 
order p 4 , good phenomenology, • • •) includes large N c QCD as well as more short- 
distance QCD constraints, see last section. 



3 Results for Bk 

Here we only recall our results in [1, 4]. We got good numerical scale matching 
for Bk for the real case and not so good for the x-hmit B K , see Figure 1. Our 
results are also analytically scheme independent [4]. 




(14) 




(15) 



B\ = -(0.43 ± 0.20) = 0.32 ± 0.15 



(16) 



in the chiral limit and 



B K = 0.77 ±0.10 
for the case of real pion and Kaon masses. 



(17) 
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Figure 1: Matching of the short-distance scale for our Bx prediction. From down 
to up curves: using CHPT at LO in the chiral limit, using ENJL in the chiral 
limit, using CHPT at LO with physical masses, using ENJL with physical masses. 
Notice the quality of this last curve. 



4 Results for e' K 

In [3] we made a prediction for e' K which we would like to update now. The new 
things we would like to input are the non-FSI corrections which after the work 
in [11, 17, 18, 19] are known. We also use the recent complete isospin breaking 
result of [6] and the new analytical results for 1m. (c 2 Ge) [5, 12, 13]. Our result 
in [3] did contain the FSI corrections but not the non-FSI which were unknown 
at that time. 

In [5, 12, 13] there are recent calculations of Im (c 2 Ge) using dispersion re- 
lations. The results there are valid to all orders in 1/N C and NLO in as- They 
are obtained using the hadronic tau data collected by ALEPH [20] and OPAL 
[21] at LEP. The agreement between them is quite good and their results can be 
summarized in 

lm{e 2 G E ) = -(4.0 ±0.9) f 8 ™ eV \ 6 Imr . ( 18 ) 
In the Standard Model 

Imr = -(6.05 ±0.50) x 10" 4 . (19) 
In [14], they used a Minimal Hadronic Approximation (MHA) to large N c 
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QCD to calculate Im (e 2 Ge) with the result 

Im(e 2 G £ ) = -(6.7±2.0) (^j'lmr, (20) 

which is also in agreement with though somewhat larger than (18). There are 
also lattice results for Im (e 2 Ge) both using domain- wall fermions [22] and Wilson 
fermions [23] . All of them made the chiral limit extrapolations, their results are 
in agreement between themselves (see comparison Table in first reference in [12]) 
and their average gives 

Im (e 2 G E ) = -(3.2 ± 0.3) ( 8 ™ eV ) 6 Imr . ( 2 l) 

The result we found in [3] for Im G 8 was 

Im G 8 = (4.4 ± 2.2) ( 8 ™ eY \ Im T ( 2 2) 
V F J 

at NLO in 1/N C . The uncertainty is dominated by the quark condensate error, 
we used [24] 

(0|gg|0) M g(2GeV) = -(0.018 ±0.004) GeV 3 (23) 

which agrees with the most recent sum rule determinations of this condensate and 
of light quark masses -see [25] for instance- and the lattice light quark masses 
world average [26]. 

We also made a calculation for Re G$ and G27 using the technique explained 
in Section 2 which reproduced the A J = 1/2 rule for Kaons within 40% through a 
very large Q 2 penguin-like contribution -see [2] for details. The results obtained 
there are 

87MeV\ 4 



ReG 8 = (6.0 ± 1.7) (— —) 



and G 27 = (0.35 ± 0.15) (^^-) , (24) 

V F ) 

in good agreement with the results obtained in [19] from a fit of K — > tttt and 
K — > 3n amplitudes at NLO to data 

/87MeV\ 4 



ReG 8 = (7.0 ± 0.6) i^— —) 
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and G 27 = (0.50 ± 0.06) ( 8 ™ eV \ . (25) 

V F J 

Very recently, using a MHA to large N c QCD, the authors of [27] found 
qualitatively similar results to those in [2, 3]. I.e., enhancement toward the 
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explanation of the AI = 1/2 rule trough large Q 2 penguin-like diagrams and a 
matrix element of the gluonic penguin Q 6 around three times the factorisable 
contribution. Indications of large values of Im G$ were also found in [28] . 

The chiral corrections to (9), (10), and (11) can be introduced through the 
Ca/=i/2, Co and C 2 factors as follows 



Re cio 


^Reao 


Re Gt2 


\Rea2 


Im ao 


( Im a 


Recto 


\Re ao 


Im a 2 


/Ima2 


Rea 2 


\Re a 2 



(2) 
(2) 
(2) 



Ca/=1/2 ; 

C ; (26) 
Im a 



C 2 + f2 c ff 



Reao 



The full isospin breaking corrections are included here through the effective pa- 
rameter f2 e ff = (0.060 ± 0.077) recently calculated in [6]. 

We get from the fit to experimental K — > irn amplitudes in [19] 

C AJ - 1/2 = # = L90±al6 = 1.22 ±0.15 (27) 
A/ " 1/2 S 2 1.56 ±0.19 V ; 

Where Si are the chiral corrections to Re aj to all orders while we call T\ to the chi- 
ral corrections to Im a/ to all orders. Therefore, they contain the FSI corrections 
which were exhaustively studied in [17] plus the non-FSI corrections which are a 
sizeable effect and of opposite direction. All these chiral corrections contain the 
large overall known factor fx/^/F^ ~ 1.47 from wave function renormalization. 

The imaginary parts Ima; get FSI corrections identical to Re ai owing to 
Watson's theorem. In addition, both due to octet dominance in Re a and Im a 
and to the numerical dominance of the non-analytic terms at NLO in K — > im 
amplitudes [11, 17, 19] 

C = %/So ~ 1.0 ±0.2, (28) 

to a good approximation. 

The situation is quite different for Ima 2 which is proportional to Im (e 2 Gs) at 
lowest order since ImG27 = in the Standard Model. From the works [17, 18, 19] 
we also know that 

T 2 0.70 ± 0.21 -O.73L4/IO- 3 
Co = — ~ 



S 2 " 1.56 ±0.18 



= 0.45 ±0.15 -0.47^^. (29) 
Putting all together, we get at LO in CHPT (using (18) and (22)) 

1 / Re a2 Im ao \ ^ 

7= 7; — -t^ — - = -(10.8 ± 5.4 Imr 

|ex|v2 VRea Rea / 
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And therefore, 



\£k 



L /^Ima,f) = (2.7±0.8)Imr. 
|V2 \Rea Re a2/ 



e K ) 



(2) 

= ((-10.8 ± 5.4) + (2.7 ± 0.8)) Imr 



(30) 



= -(8.1 ± 5.5) Imr = (4.9 ± 3.3) x 10~ 3 . (31) 

This result is scheme independent and very stable against the short-distance scale 
as can be seen in Figure 2. The difference with the result in Figure 2 [3] is due 




0.004 



0.002 



1 1 1 1 1 1 

0.5 0.6 0.7 0.8 0.9 1 

H [GeV] 

Figure 2: Matching of the short-distance scale dependence of our LO in CHPT 
e' k /sk prediction. Labels I and II are for two different values of as- The two 
curves for two choices of perturbative matching, see [3]. Notice the quality of the 
matching. 



to the new values of Imr in (19) and Re(?8 and G27 in (25). 

Including the known and estimated higher order CHPT corrections, we get 

1 Rect2 Ima 



\e K \y/2 Rea Rea 

and 

1 Rea 2 Ima2 



-(8.9 ±4.8) Imr (32) 



I £ a- I -\/2 Rea Rea 2 
((1.0 ± 0.3) + (0.5 ± 0.7)) Imr = (1.5 ± 0.8) Imr . 



(33) 
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where the second part comes from isospin breaking contribution with Q e ff = 
0.06 ± 0.08 [6], and we used L 4 = in (29). And therefore, 




= ((-8.9 ± 4.8) + (1.5 ± 0.8)) Imr 
(7.4 ± 4.9) Imr = (4.5 ± 3.0) x 10" 3 



(34) 




(1.66 ±0.16) x 10~ 3 . 



(35) 



Though the central value of our Standard Model prediction in (34) is a factor 
around 3 too large, within the big uncertainties it is still compatible with the 
experimental result. Two immediate consequences of the analysis above, namely, 
the LO CHPT prediction (31) is actually very close of the the final result (34) 
and second, the part with AI = 1/2 dominates when all higher order CHPT 
corrections are included. 

5 Conclusions and Prospects 

The large final uncertainty we quote in (34) is mainly due to the uncertainties 
of (1) the chiral limit quark condensate, which is not smaller than 20%, (2) L 5 , 
which is around 30%, and (3) the NLO in 1/N C corrections to the matrix element 
of Q e , which is around 20%. All of them together make the present prediction 
for the AI = 1/2 contribution to e' K to have an error around 55%. Reduction in 
the uncertainty of all these inputs, especially of the quark condensate and L5 is 
needed to obtain a reasonable final uncertainty. 

We substituted the value used in [3] for lm(e 2 G E ) by the one in (18), notice 
however that numerically they coincide within errors. There are also large uncer- 
tainties in the AI = 3/2 component coming from isospin breaking [6], and more 
moderate in the non-FSI corrections to Ima2, fortunately the impact of them in 
the final result is not as large as the ones associated to the AI = 1/2 component. 

Assuming the AI = 3/2 component of e' K is fixed to be the one in (18) as 
indicated by the analytic methods [5, 12, 13, 14] and lattice [22, 23], one can try 
to extract the value of ImG 8 from the experimental result in (35). We get 



i.e. the central value coincides with the large N c result for ImGg within 30% of 
uncertainty. This value has to be compared to our result in (22). 

We have presented in [31] a ladder resummation hadronic model which keeps 
the good features of the ENJL model -CHPT at NLO, for instance- and improves 




(36) 
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adding more short-distance QCD constraints and the analytic structure of large 
N c . We intend to use the Green's functions calculated within this hadronic model 
analytically, and study the relevant internal cancellations and dominant hadronic 
parameters in the quantities we calculate with them. Work in this direction where 
we will study the origin of the large chiral corrections to Bk we found in [1] is in 
progress [32]. This program will also be extended to study the origin of the large 
value for ImG 8 in (22) which we got in [3]. Large values of ImG 8 were previously 
pointed out in [28] and more recently in [27]. 
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